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© A system for predicting the life of a cutting tool in an automated metal cutting machine, comprises a plurality 
of sensors for measuring physical data (601) relating to a cutting operation of said machine, and a computer 
connected to these sensors. The computer comprises a feature extractor (603) for transforming the physical data 
into feature values, and a tool life predictor (605) for arriving at a prediction of the life of the cutting tool. The tool 
life predictor has an influence diagram having an input level of feature nodes, intermediate nodes connected to 
the feature nodes and to a machining class node, and an output node connected to the intermediate nodes. The 
tool life predicator further comprises a first module for computing averages of feature values for each sensor, a 
second module responsive to the feature averages for making a classification of the feature average as indicative 
of a long life or a short life, and a third module connected to the influence diagram for combining the 
classifications to arrive at a classification and a conditional probability at each intermediate node and a 
classification at the output node which is the predicted life of the cutting tool. 

Furthermore, a method to enable a real-time expert system to learn is provided which comprises the steps of 
creating a relational structure corresponding to the influences that a first and a second variable have on a third 
variable, obtaining data for these variables over some duration of time, determining a tolerance zone around a 
sample which divides classes and using a dynamic distance classifier, which is responsive to the tolerance zone, 
to determine decision boundaries for the first and second variables and probabilities associated with the decision 
boundaries. The probabilities indicate the confidence with which a decision boundary classifies values of the first 
and second variables with respect to the value of the third variable. 



Q- 
LU 



Rank Xerox (UK) Business Services 




EP 0 509 817 A1 



START DRILLING 
AND 
MEASUREMENTS 



601 



603 



COMPUTE FEATURES FROM 
FIRST THREE HOLES OF OATA 



605 



PREDICT LIFE OF ORILl 
rt LONG LIFE OR SHORT LIFE 




CONTINUE DRILLING 
AND 

VERIFY AT MID-LIFE 




6/3 



OBTAIN 
REPLACEMENT 



6/7L 



CONTINUE DRILLING 
UNTIL PREDICTED 
END OF LIFE 



6/9 



INITIATE 
CONTINOUS 
OIAGNOSEO 



62L 



MEASURE AND COMPUTE 
FEATVRES FROM NEXT 
THREE HOLES 



-l 625 



CONTINUE 
(WILLING 

UNTIL 
PREDICTEO 
END OF LIFE 



CONTINUE 
DRILLING UNTIL 
TOOL DIAGNOSED WORN 



I 



REPLACE 
TOOL 



623 



F/o.3 



2 




EP 0 509 817 A1 



NOTICE 

(C) Copyright Texas Instruments Incorporated 1991. A portion of the disclosure of this patent document 
contains material which is subject to copyright protection. The copyright owner has no objection to the 
5 facsimile reproduction by anyone of the patent document or the patent disclosure, as it appears in the 
Patent and Trademark Office patent file or records, but otherwise reserves all copyrights whatsoever. 

BACKGROUND OF THE INVENTION 

10 1 . Held of the Invention 

This invention relates to manufacturing automation and in particular to a system and method for timely 
replacement of cutting tools through early prediction of cutting tool life and accurate diagnosis of tool wear 
using an Influence Diagram Expert System and multiple sensors. 

75 

2. Description of the Related Art 

The background of this invention is described in conjunction with manufacturing automation as an 
example. 

20 Automation in manufacturing is a rapidly growing field with significant impact on improved product 
quality, high productivity and reduced manufacturing cost. While it has been adopted successfully at 
different levels of manufacturing process, one of the main hindrances of completely automating any metal 
cutting operation is the timely replacement of cutting tools. Currently, human intervention is required to 
change the cutting tools at appropriate heuristically chosen intervals and typically, one machinist is 

25 assigned to supervise and replace cutting tools for a batch of automated machining centers. With adequate 
sensors for in-process monitoring of the tool condition, an improvement can be obtained in machining 
economics, especially in high volume production lines. Cutting tool condition monitoring, wear diagnostics 
and appropriate machine control problems have been investigated through different methodologies: real 
time expert system based approaches and mathematical model based approaches. While these on-line 

30 monitoring techniques provide a solution to correctly identify a damaged or worn tool, they are postmortem 
techniques; waiting for an appropriate decision after the tool is worn or damaged. 

A tool changing operation may involve multiple tasks like procurement of a replacement tool from a 
centralized storage area, loading tool-setting software from the plant manager to the machine controller, 
dynamic re-scheduling of the production line to account for the anticipated machine down time, and so on. 

35 The present tool monitoring techniques do not provide adequate lead-time to initiate these multiple tasks in 
preparation for a tool-change, and hence, no significant improvement is made in reducing the downtime due 
to cutting tool changes. Another added limitation of present continuous on-line monitoring is the requirement 
of dedicated processors on each machine, primarily because of stringent response time requirements of the 
present day real time diagnostic systems. 

40 In order to automate the cutting tool replacement task with adequate forewarning, a reasonably accurate 
estimate of the life of the cutting tool is required. As mentioned earlier, a significant advantage would be 
that the predicted tool life information could be utilized for dynamically re-scheduling machining operations 
and cutting tool replacements with a reduction in inventory and labor costs. 

Mathematical model based techniques, called Taylors tool life equation, and extensions thereto, may 

45 be applied to estimate the cutting tool life. However, these equations provide a poor estimate of the life of 
the cutting tool. For example, in the case of drilling, work piece material hardness is an important factor 
affecting the life of the drill but is not taken into account. While modified forms of Taylor's equation with 
corrections for work piece hardness have been applied to determine the drill tool life, such an approach is 
impractical in a production line since hardness for each work piece is not readily available or easily 

so measurable. Furthermore, in some other machining situations, Taylor's model has been found unsuitable. 

Therefore, it would be desirable to have a system in which a cutting tool's life can be accurately 
predicted during the initial use of the tool so as to provide adequate lead-time for the replacement of the 
tool. 

In the area of cutting tool wear diagnosis, some present day techniques are pattern recognition (PR), 
55 neural nets (NN) and real-time expert systems (RTES). PR and NN, use pattern samples that are already 
classified based on some wear criteria and train the system to recognize this fact. This approach is called 
'supervised learning'. In 'unsupervised learning', no wear criteria is used and hence, the knowledge about 
the class to which the samples belong is not available a priori. Such an unsupervised training approach 
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using fuzzy clustering technique has been applied to the drill tool classification problem. Using thrust and 
torque data from a torque dynamometer, a fuzzy classification of the state of the drill tool has been 
performed successfully. However, the drill is classified as worn only after the cutting edges are severely 
damaged. This prediction is delayed until the drill is about to fail. 

5 The advantages of the RTES techniques are the ease of modification, the ability to select optimal 
control decisions by optimizing over a cost function after classifying the current state of the tool, and the 
facility to repetitively utilize the same inference engine while applying to different machining operations. If 
the diagnostic system needs to be updated to handle new combinations of machining parameters or new 
sensors, then in case of PR or NN techniques, the system has to be retrained with new as well as old data, 

70 a rather time consuming process. However, in case of RTES, the knowledge base can be easily updated 
with the new information. 

Heretofore, one of the problems with the present RTES technique has been the lack of learning 
capability. The development and tuning of the RTES requires detailed analysis of pertinent data to extract 
features and subjective estimate of the conditional probability distributions that form the critical elements of 
75 the knowledge base. 

SUMMARY OF THE INVENTION 

A system for predicting the life of a cutting tool in an automated metal cutting machine, comprises a 

20 plurality of sensors for measuring physical data relating to a cutting operation of said machine, and a 
computer connected to these sensors. The computer comprises a feature extractor for transforming the 
physical data into feature values, a tool life predictor for arriving at a prediction of the life of the cutting tool. 
The tool life predictor has an influence diagram having an input level of feature nodes, intermediate nodes 
connected to the feature nodes and to a machining class node, an output node connected to the 

25 intermediate nodes, a first module for computing averages of feature values for each sensor, a second 
module responsive to the feature averages for making a classification of the feature average as indicative of 
a long life or a short life, and a third module connected to the influence diagram for combining the 
classifications to arrive at a classification and a conditional probability at each intermediate node and a 
classification at the output node which is the predicted life of the cutting tool. 

30 Furthermore, a method to enable a real-time expert system to learn comprises the steps of creating a 
relational structure corresponding to the influences that a first and a second variable have on a third 
variable, obtaining data for these variables over some duration of time, determining a tolerance zone around 
a sample which divides classes and using a dynamic distance classifier, which is responsive to the 
tolerance zone, to determine decision boundaries for the first and second variables and probabilities 

35 associated with the decision boundaries. The probabilities indicating the confidence with which a decision 
boundary classifies values of the first and second variables with respect to the value of the third variable. 

BRIEF DESCRIPTION OF THE DRAWINGS 

40 In the drawings: 

FIG 1. is a block diagram of a computer aided manufacturing (CAM) system; 

FIG 2. is a schematic diagram showing the logical links between sensor data in a CAM system and an 
Influence Diagram Expert System (IDES). 

FIG 3. is a flow chart for operating a CAM system using an IDES to predict tool life and to monitor tool 
45 wear; 

FIG 4. shows the relationship of the training phase to the operations phase; 
FIG 5. shows an example influence diagram; 

FIG 6. shows the influence diagram topology for predicting tool life according to the present invention; 
FIG 7. is a graph of experimental data used during the training phase of the present invention; 
so FIG 8. shows the influence diagram topology for monitoring tool wear according to the present invention; 
FIG 9. is a flow chart showing how a dynamic distance classifier according to the present invention is 
utilized in a computerized machining operation. 

Corresponding numerals and symbols in the different figures refer to corresponding parts unless 
otherwise indicated. 

55 

DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS 

In Figure 1, is shown in an automated manufacturing operation using this invention. A metal cutting 
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machine 90 has a drill bit 100 attached to a quill 108. rotated by a spindle motor 106. The spindle motor 
106 is connected to an AC electrical power source 117 via wires 116. 

Drill bit 100 is lowered by feed motor 110 into a work piece 102 attached to a work table 126 by vise 
104. Feed motor 110 is connected to a DC electrical power source 119 via wires 120. 

5 A number of sensors measure physical data pertinent to the drilling operation of machine 90. A spindle 
motor current gage 114 senses the AC current transmitted through wires 116, a feed motor current gage 
118 measures the DC cuttent transmitted through wires 120, a strain gage 122 is attached to quill 108 for 
measuring the thrust force exerted by the tool 100 on work piece 102, an accelerometer 128 senses 
vibrations in vise 104, and a dynamometer 124 measures the force exerted by vise 104 on table 126. The 

io outputs from these sensors are fed into an analog signal conditioner 130 for adjusting the raw signals into a 
range acceptable as input by analog-to-digital converter 132. Conditioned signals from the analog signal 
conditioner 130 are transmitted through lines 160 through 170 into analog-to-digital converter 132 of a 
digital computer 134. Computer 134 contains a central processing unit 136 and a memory 138. Connected 
to computer 134 are monitor 144, keyboard 148, disk drive 140, and tape drive 142. 

75 Adjacent to machine 90 is a measurescope 148 which has a monitor 150 attached to it. The 
measurescope 148 is positioned so that it can visually measure the wear scars of drill bit 100. 

Computer 134 processes inputs from sensors 114, 118, 122, 124, and 128 to predict when drill bit 100 
will be worn out and at the end of predicted life to monitor drill bit 1 00 for wear on a continuous basis. The 
prediction and diagnosis of wear may either be displayed on monitor 144, for a human operator to take 

20 some action, such as replacing the tool, or may initiate the automatic procurement and replacement of drill 
bit 100. The replacement may be done by some process control hardware 112 attached to computer 134. 

During certain phases of the operation of machine 90, computer 134 will be idle with respect to 
processing data and controlling machine 90. That idle time may be used to service other machines (not 
shown). Thus, computer 134 may be used to predict tool life and diagnose tool wear in multiple machines 

25 which operate in parallel. 

The operation of machine 90 is controlled by an Influence Diagram Expert System (IDES) 154 shown in 
FIG 2. As Figure 2 shows, computer 134 executes Influence Diagram Expert System (IDES) 154. It predicts 
the length of life of drill bit 100, monitors drill bit 130 for wear, and initiates the replacement procedure 
when drill bit 100 has been diagnosed as worn out. Figure 2 shows a schematic of the logical links between 

30 sensor data 156 and IDES 154. The sensor data is, for example, sensors 114, 118, 122, 124, and 128, of 
Figure 1. IDES 154 is also connected to knowledge base 152, in which IDES 154 stores information it uses 
to process data obtained from the sensors to make decisions which are displayed on monitor 144 or which 
initiate some action on the part of process control hardware 112. Knowledge base 152 may be physically 
located in memory 138 of Fig. 1, on disk 140, on tape 142, or distributed over these devices or other 

35 storage devices not shown. 

IDES 154 is also logically connected to parameters 156. These may, for example, be entered from 
keyboard 148, of Fig. 1, or retrieved from data files on disk 140 or tape 142. Parameters 156 includes, for 
example, information for categorizing the machining operation, such as, tool size, work piece material, drill 
feed rate, and drill speed. 

40 The operation of machine 90 and IDES 154 consists of the following phases: 

1. An experimental phase in which the drilling operation is carried out and sensor measurements are 
taken at regular intervals throughout the life of the tool. Furthermore, the point at which the drill has been 
worn out, i.e., its life has been expended, is noted. As each hole is drilled its number is noted and each 
data point is indexed by it. Similarly, the point at which the tool is deemed to have been worn out is also 

45 noted by hole number. This experimental data is later used to build knowledge base 152 for both tool life 
prediction and tool wear diagnosis. The knowledge base 152 may be stored on disk 140 or tape 142. 

Both raw sensor data and quantities computed from the raw sensor data are used in both tool life 
prediction and tool wear diagnosis. The manner in which the computed quantities are derived from the 
raw sensor data is discussed in greater detail below. The raw sensor data and the computed quantities 

50 are collectively called features. 

2. A training phase, in which the knowledge base is built and fine-tuned consists of two distinct aspects: 
first, building the knowledge base for accurate prediction of tool life and, second, the knowledge base for 
the diagnosis of tool wear. Certain portions of the training phase are automated by using the Dynamic 
Distance Classifier algorithm described below. 

55 3. An operations phase in which IDES 154, utilizing the knowledge base which was built during the 
training phase, predicts tool life and monitors tool wear. 

Thus, the goal is to build a knowledge base which an IDES 154 uses to predict tool life and to monitor 
drill wear. The training phase is more easily understood if the operations phase is considered first. 
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OPERATIONS PHASE 

FIG 3. is a flow chart showing the operation of a computer aided manufacturing operation according to 
the present invention. A drilling operation starts at Step 601 . In conjunction with the drilling, measurements 

5 are taken from the various sensors used in the Influence Diagram to predict tool life. Early in the operation, 
for example after three holes, the features are computed, Step 603. These features are used by IDES 154 
to predict the life of drill 100, Step 605. 

If the prediction is that the tool will have a long life (Step 607), the machine continues drilling until the 
tool's predicted mid-life, Step 609. At which point, a verification procedure is invoked. If the earlier 

io prediction of a long life is verified the drilling continues until the end of the predicted life of the tool, Step 
615. At that point sensor data is collected and features computed for the next three holes, Step 619. The 
process cycles back to predicting the remaining life of the drill, Step 605. 

During the time between prediction of a "long life" (Steps 605 and 607) and the verification at mid-life 
(Step 609) computer 134 does not process data from machine 90. Thus, during that time, computer 134 is 

is available for other tasks, for example, servicing other machines attached to it. 

If a long life is not predicted in Step 605, or if a predicted long life is not verified at the predicted mid- 
life in Step 609, then a replacement tool is obtained, Step 613. One embodiment is represented by the 
dashed path from Step 613 to Step 623. Drilling continues until the predicted end of the life of the tool. Step 
625. The tool is replaced at the point at which IDES 154 has predicted that the tool's life ends, Step 625. 

20 An alternative embodiment is indicated by the solid path emanating from Step 613. Because it is 
desirable to never use a tool beyond the point at which it is worn out, the tool life prediction is adjusted to 
be a conservative estimate. For that reason, tools which have not worn out may have been predicted as 
worn out. Rather than discarding such tools, the embodiment represented by the path through Steps 613, 
617, and 619, provides a method for using the sensors 114, 118, 122, 126, and 128, to diagnose the tool for 

25 wear. Thus, IDES 154 continuously monitors the sensors and diagnoses tool wear, Step 617, and drilling 
continues, Step 621, until IDES 154 diagnoses the end of the life of the tool, Step 621. When the tool has 
been expended the tool is replaced, Step 623. The entire operation starts over with the fresh drill at Step 
601. 

. The method described above in conjunction with Fig. 3 may be implemented to be executed on 
30 computer system 134. One embodiment has a functional module to compute features from the raw sensor 
data (Step 603), a functional module each for the tool life prediction (Step 605), the verification of the 
prediction (Step 609), and the tool wear diagnostics (Steps 617 and 621). 

TRAINING PHASE 

35 

Fig. 4 shows the relationship of the training phase for tool life prediction to the operation of the IDES 
154 and the automated manufacturing system shown in Fig. 1. A similar approach is used for training the 
system for tool wear diagnosis. 

The training operation may be divided functionally into two major categories: off-line knowledge base 

40 generation and on-line tool-life prediction. The off-line task generates the qualitative and the quantitative 
knowledge base for the influence diagram model of the drilling process. The knowledge base is used by the 
real-time system for on-line prediction of tool-life and tool wear diagnostics. The input to the system are the 
multi-sensor raw data and the cutting tool wear measurements collected during the experimental phase as 
well as the machining category. The term "machining category" represents the combinations of machining 

45 parameters that are used for collecting the training data. When the machining parameters vary considerably, 
the sensor features also vary widely. Hence, without any knowledge about the machining conditions and 
with just the sensory data it is difficult to establish a diagnostic mapping between the sensory features and 
the cutting tool state. Subsequently, the knowledge engineer divides the broad spectrum of machining 
conditions into several categories such that each category represents the machining parameters that 

50 resemble close to each other. An example of machining categories is given in Table 1. 
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Category 


Workpiece 
material 


Orill 
size/type 


Speed 
(rpro) 


Feedrate | 
(inch/rev.) f 


0 


Ductile cast 
iron 


0.25* 
Tungsten 


1100-1400 


0.006-0.011 


1 


Ductile cast 
i ron 


0.25" 

Tungsten 


1600-2000 


O.C03-0.006 


2 


A36 steel 


0.25" High 
Speed Steel 


3600-3900 


0.002-0.004 


3 


A36 steel 


0.25* High 
Speed Steel 


2300-2700 


0.002-0.004 


4 


A36 steel 


0.25" High 
Speed Steel 


2300-2700 


0.008-0.010 


5 


A36 steel 


0.125" High 
Speed Steel 


3700-4000 


0.006-0.007 


6 


A36 steel 


0.125' High 
Speed Steel 


3700-4000 


0.009-0.010 


7 


A36 steel 


0.125" High 
Speed Steel 


2900-3100 


0.009-0.010 



Table 1. Categorization of Widely Varying Machining Conditions 

The knowledge base generation is divided into several sub-tasks: feature selection, building the 
30 relational structure, determining feature state descriptors, and subjective estimation of the conditional 
probability distributions. 

1 . Feature Selection 

35 In Fig. 1, the spindle motor current sensor 114 is an AC current sensor. It is a noninvasive induction 
sensor that is simply clipped around the input power line. The AC current sensor output is rectified and low 
pass filtered in order to acquire low frequency energy of the sensor signal. This energy is directly 
proportional to the cutting torque exerted by the tool on the work piece. As the tool wears, the cutting force 
experienced by the tool increases. This increase in cutting force causes a corresponding increase in spindle 

40 motor current. The RMS value of spindle motor current thus becomes a valuable feature for observing the 
trend in tool wear during the initial phase of drilling. 

In addition to the RMS value, the change in RMS value compared to a fresh tool (ARMS) is a useful 
feature as it indicates the temporal trend of the cutting torque. A faster rate of change indicates a rapid wear 
and hence shorter tool life. 

45 Another valuable source of information about the drilling process is the spindle mounted strain gage 
122. It indicates the thrust force exerted by the tool 100 on the work piece 102. As the tool 100 wears, the 
thrust forces increase due to a constant feed rate and an increase in energy required to penetrate the work 
piece 102. The actual value used is the mean of the strain gage force during the drilling of one hole. 

Furthermore, as in the case of the spindle motor current sensor 1 1 4, the raw value (the mean) as well 

so as its incremental change from a fresh tool (Amean) are used to predict the tool life. 

In the preferred embodiment the features extracted from the measurements of the spindle motor current 
sensor 114 and the strain gage sensor 122 are used to predict tool life and to continuously diagnose tool 
wear at the end of the predicted life of the tool. However, in an alternative embodiment other sensors, for 
example the feed motor current sensor 110, the accelerometer 128, and the dynamometer 124, are used 

55 either in lieu of or in conjunction with either the spindle motor current sensor 114 or the strain gage sensor 
122. 

2. Building the relational structure 
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An Influence Diagram Expert System (IDES) is an expert system in which relationships between 
variables are qualitatively and quantitatively represented and manipulated to arrive at some decision. An 
influence diagram, i.e., the relational structure, which is a graphical network and theoretical structure for 
knowledge, provides probabilistic inference and expected value decisions, and has two hierarchical levels: 

5 symbolic and quantitative. At the symbolic level, an influence diagram consists of an acyclic directed 
network with nodes representing relevant system state variables and the arcs representing state variable 
interrelationships. For example, the diagram in Fig. 5 is a simple influence diagram which schematically 
shows the probabilistic influence three variables, RMS, ARMS, and category, have on the state of a cutting 
tool (wearl). At the quantitative level the mapping between the relevant states are represented by discrete 

10 conditional probability density functions. For example, in the sample influence diagram in Rg. 5, node 3 
would be quantitatively represented by the conditional probability { wearl \ RMS, ARMS, category }. Also, 
in a general model, the root nodes 1 , 2, and 3 are quantitatively represented by their respective marginal 
probabilities. However, when the influence diagram technique is used in a real-time paradigm, they are 
represented by uniform distributions. A diagnostic inference, for example, the marginal probability of the 

15 cutting tool being worn, results from a set of transformations applied to the diagram in response to a 
specific query. The topological structure remaining after the set of transformations represents the answer to 
the query in symbolic form. At the quantitative level, these transformations would provide the answer to the 
query in numerical form. Three types of transformations used in this reduction process are Sensor Node 
Removal (table look-up operation), State Node Removal (applying summation rule of probability) and Arc 

20 Reversal (applying Bayes' rule). Once diagnostic inference has evaluated the likelihood of critical events, 
the best course of action must be determined. 

Fig. 6 is an influence diagram for predicting tool life in the preferred embodiment. Verification uses an 
influence diagram with the same topology. However, at the numerical level, the state space of each feature 
and the subjective estimates of the probability distributions are different for different phases. Hence, the 

25 portion of the knowledge base corresponding to the influence diagram for the verification phase is different 
from the corresponding portion of the prediction phase. 

The features Arms and rms from the spindle motor current sensor influence the //Ye/ variable, an 
intermediate node in the influence diagram, strain gage force features Amean and mean influence the Iife2 
variable, and machining parameters speed, feed, work-piece material and tool size influence the class 

30 variable. The latter is not a probabilistic influence but rather a categorization. An example categorization 
based on these particular parameters is shown in Table 1 . However, other categorizations are possible. The 
class variable also influences the lifel and Iife2 variables. The range of values for life 1 and Iife2 is "long 
life" or "short life." Thus, given particular values of each feature, for example rms and Arms, both values 
together indicate with a particular conditional probability that the value of life 1 is "long life." 

35 Finally, life 1 and Iife2 influence the life variable, which is an output node in the influence diagram. The 
values available for the life variable is "long life" and "short life". Both of which are values dependent on 
the type of tool in question. For example, for 0.25" drills, based on experimental results, "long life" is 
defined as 100 holes and "short life" as 10 holes, and similarly for 0.125" drills, "long life" corresponds to 
50 holes and "short life" to 7 holes. However, in alternative embodiments these definitions may be different. 

40 

3. Determining feature state descriptors 

The quantitative level of the influence diagram, as part of knowledge base 152, each node is divided 
into many states depending on the range of values, symbolic or numerical, the variable can take. The 

45 values that divide the range into many states are referred to as the state descriptors. The number of states 
varies. For example, for the fife node 423, the lifel node 417, and the Iife2 node 419, there are two states 
"long life" and "short life." The class node 421, has as many states as the total number of categories, e.g. 
eight in the case of Table 1. In one embodiment, using the categories in Table 1, node 415 has "0.25-inch" 
and "0.125-inch" states, and node 413 has the states "steel" and "cast iron." Thus, the states for nodes 

50 423, 417, 419, 415 and 413, are all qualitative. The remaining nodes in Fig. 6 have state descriptors defined 
by quantitative intervals. The states for nodes 409 and 41 1 are obtained by dividing the normal operating 
ranges of speeds and feeds respectively, over as many intervals as needed to uniquely define each 
category in Table 1 . For node 41 1 , and Table 1 , a state 0 represents feed rate less than 0.008 ipr, and a 
state 1 represents feed rates greater than or equal to 0.008 ipr. Similarly, in the case of node 409, a state 0 

55 is defined by drill speeds less than 3200 rpm, and a state 1 is defined as speeds greater than or equal to 
3200 rpm. In the case of feature nodes 401 through 407 the state descriptors are decided by analyzing and 
correlating the raw experimental data with the actual life for the tool for all the categories. 

Fig 7. is an illustration of raw data collected for two category 5 cases. It illustrates how to obtain the 
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decision boundaries associated with the RMS feature. Thus, the figure represents a specific feature, RMS, 
of node 403, and a specific category, category 5. The experimental data may be plotted in similar graphs 
for other features and categories. In Fig. 7, the RMS feature of spindle current signal for both drills is plotted 
as a function of number of holes drilled. Each point in the graph represents the RMS value for signals 

5 collected for one hole. The vertical dotted lines (labeled "Drill #1 Worn" and "Drill #2 Worn", respectively) 
indicate the hole number at which the drills were actually declared as worn by wear measurement criterion. 
This hole will be referred to as the "wear hole." The circles to the left of vertical line for drill #1 belong to 
the good tool condition and ones to the right belong to the worn tool condition for drill #1. Similarly drill #2 
is represented by plus symbols. Note that the life of Drill #1 is a "long life" whereas the life of Drill #2 is a 

70 "short life." The objective is to select one threshold (horizontal line called decision boundary) value that 
separates the "long life" state from the "short life" state for both the drills. For the prediction phase the 
average RMS value for the first three holes is used. Thus, the decision boundary for the prediction phase 
(labeled a) is selected so that the average RMS value for the first three holes will indicate whether a "long 
life" is indicated or not. In the example of Fig. 7 the average RMS for the first three holes of Drill #1 is 

75 approximately 0.45 and the average RMS for the first three holes of Drill #2 is approximately 0.36, therefore 
the selected decision boundary for RMS for Category #5 drills is 0.4. Similar analysis is done for the 
verification phase. However, in the case of verification the maximum feature values for the drills up to the 
mid-life is used rather than the average for the first three holes. In the case of Drill #1 the max value is 0.59, 
as indicated by the boxed ring. The max value for Drill #2 is 0.5. Therefore, the selected decision boundary 

20 for the verification phase for RMS for category #5 drills is 0.57. 

The decision boundaries for all features and all categories is done in the same manner as explained 
above for RMS and category #5. Thus, a quantitative knowledge base of decision boundaries is built. 



4. Subjective estimation of the conditional probability distributions. 

25 

The next step is to assign marginal probabilities associated with the decision boundaries. The marginal 
probability for a given feature and category is a subjective evaluation on the part of the knowledge engineer. 
It indicates the probability that the decision boundary accurately predicts "long life." 

Having assigned the marginal probabilities for all feature nodes 401 through 407, the knowledge 
30 engineer determines conditional probabilities for the intermediate nodes 41 7 and 41 9. The class node, node 
421, representing machining category, is assigned a binary deterministic conditional distribution (0 or 1) 
given any element of its conditional cross product space consisting of all the combinations of states of the 
predecessor nodes 409 through 415. 

Using node 417 as an example, for each of the predecessor features rms (node 403) and Arms (node 
35 401) a subjective weight on a scale of 0 to 1 is assigned based on how well each decision boundary 
separated the "short life" from the "long life" states of all the categories. Let these weights be tV403 and 
tv*oi for rms and Arms, respectively. For each category, critical states are identified for nodes 401 and 403, 
below which the feature values correspond to the "long life" state and above which the feature values 
correspond to the "short life" state. For some category c, let $ao3 and s*oi, be the critical states for nodes 
40 403 and 401 respectively. Let 1Uo3 and O401 be the state spaces for nodes 403 and 401 respectively. The 
conditional probability assessment is made as follows. 

For category c, the conditional probability, P, that the node 417 is in "short life" for (/ x j) c 0*03 x CUoi 

P=max(w 403 ,w 401 ) +A t (Q-i) ; if i*s 403 , j 2S 401 
P* ^ — +A 2 (Q) ; if izs 4 o2<3 <s 4oi 

403 w 40l 

P- — — — 3 (Q) ; if i <S4O3>J* 3 40t 

P=O.Q; if i<s 40J ,j<s 401 



45 



50 



where 

55 

Q = {i + j + S403 - s*oi + 1) 



and Ai, A 2 , and A 3 are smoothing factors that gradually change the conditional probabilities in discrete 
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steps as the features' states gradually change from those corresponding to a drill with a short life to those 
corresponding to a drill with a long life. They typically assume values ranging from 0.05 to 0.1. and are 
chosen subjectively. All these subjective quantities are different for different categories and this procedure 
is repeated for all the categories. While tuning the knowledge base, they are iteratively updated until the 

5 system responds correctly for the training data. Typically no more than five iterations are necessary. The 
equations above are heuristic, and a person skilled in the art will realize other possible equations. 

As to the training phase with respect to tool wear diagnosis, Fig. 8 shows the topology of the influence 
diagram used by IDES 154 during the continuous diagnosis of drill wear in Step 617, of Fig 3. A comparison 
with Fig. 6 indicates the similarity between the influence diagram for tool life prediction and tool wear 

70 diagnosis that takes place at the end of the predicted life of the tool. However, interior and output nodes 
differ. In the tool wear diagnosis influence diagram, the features derived from the spindle motor current 
influence the wear/ variable and those derived from the strain gage force influence the weBr2 variable. 
Those intermediate nodes influence the output node, wear. The values for wear/, wear2 and wear are 
"Fresh drill" and "Worn drill." Furthermore, at the numerical level, the decision boundaries and probabilities 

75 differ from the influence diagrams associated with tool life prediction. 

In the case of tool life prediction discussed above the knowledge base is manually built by a knowledge 
engineer. The knowledge base for the continuous diagnosis of tool wear (of Step 621 of Fig. 3) is 
automatically created using a dynamic distance classifier. Fig, 9 shows the relationship of the dynamic 
distance classifier to the total operation of the computer aided manufacturing system having an influence 

20 diagram expert system. 

The dynamic distance classifier generates the decision boundaries between the fresh and worn tool 
classes as well as the probability of correct classification for each feature. It uses the Euclidian distance 
between the sample points and the mean of each class as the measure for classification. It bears the name 
"dynamic" due to an iterative update of the fresh tool class mean based on a new separability criterion 

25 defined below. The method is generic in that it is applicable to drilling as well as other machining 
operations, for example milling or turning. Furthermore, applications outside of the computer aided 
manufacturing environment will be apparent to a person skilled in the art. 

The method assumes that the data is collected at regular intervals for a short period of time and the 
sample points, classified based on some tool wear criteria, represent a time history of the state of the tool. 

30 This assumption is not a restriction but is in-line with the true situation in a machining environment. 

1 . Nomenclature 



The variables used in the method are defined as follows: 



35 


Co- 


Fresh cutting tool class 




Ci - 


Worn cutting tool class 




d- 


Decision boundary 




w - 


Actual sample # that divides classes 




T- 


Tolerance zone, in terms of number of samples, about w (e.g., T = {4,5,6,7} given w = 5) 


40 


W m ln " 


Minimum of / for all /"« T (e.g., w min = 4 given T = {4,5,6,7}) 




X - 


Feature under evaluation for its classification capability 




*\ - 


Value of x at / 1h sample 




9i- 


The feature mean of class Q 




P- 


Probability of successful classification for the feature under evaluation 



45 

2. Separability Condition for Cutting Tool Classes 



Conventionally, in pattern classification using linear discriminant function F, two classes Co and G are 
separable if 
50 F(xO < 0 for all X\ e Co and 
F(x f ) > 0 for all *i t Ci , 

where X| denotes the feature value for sample / and the decision boundary d is defined by 
F(d) = 0. 

Since it is undesirable to have any fresh tool states to be classified as worn, the condition 
55 F(x,) < 0 for all x t t Co 

is a requirement for correct classification. After w samples the tool is defined as worn based on some wear 
criterion. For conventional separability, the classifier diagnoses the tool to be worn for each sample 
collected beyond the w* h sample. However, in case of robust diagnosis, at the first instance of wear 
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diagnosis, a decision would have been made; either to replace the cutting tool or continue operation in 
degraded mode, depending on the economics of that particular machining operation. Hence it is sufficient 
that there exists just one correct classification of the worn drill within a tolerance zone (T) of a few samples 
about the actual worn tool sample w. The tolerance zone T may consist of a few fresh tool samples before 
5 the worn tool sample w (as indicated in the definition of 7). Consequently, fresh and worn tool classes are 
defined as separable if 

1 . F(xO < 0 for all x t € Co and 

2. there exists i e T such that F(xQ > 0, i.e., x, * G 

This relaxed but adequate separability condition simplifies the search for a decision boundary that 
to would provide robust classification. For each individual feature, this decision boundary is nothing but a 
point. This one dimensional search for the decision boundary is considerably faster than the search for the 
n-dimensional decision surface where the features are arranged into an n-dimensional vector. 

3. Classifier Method 

15 

The following pseudocode determines the decision boundary and the marginal probability for each 
feature. 

For (x = Feature(l) to Feature(n) ) 

20 

Begin For loop 1 

(1) Compute mean of each class (g } ; j = 0,1). 
25 (2) If (g 0 > g x ) then return (p = 0) and exit. 



30 
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(3) For (i = First sample to Last sample ) 
Do classification using decision rale, 

x t € C k if |x, - g k \ - min |x, - j - 0,1 

(4) If (x, e C 0 for all i e T) then return (p = 0) and exit, 

(5) If (x x € C Q for all i < w m .J and 

(there exists i e T such that x i € C x ) then 

return d = 

2 

and exit. 

(6) While ( there exist i € T such that x, € C { ) 

and (there exists i e T such that x, e C x ) 
BEGIN While 



1 - 



2 

S 



v - 1 

25 aln * 



9 0 - 



where s - \S\ , 

S - {xj ! € Ci for aiJ i < w^} 



S ♦ 1 

where Xj e C x for aii i < 



If (g 0 * g x ) then exit While loop. 
Repeat step (3) only. 
END While 

(7) If (x s e C 0 for all i < w ui J and (there exists i € T 
such that x, € C x then 

45 returned = 9 ° * 9l and p = 1) ajid exit. 



else return (d,p) and exit. 
END For loop 1 



Thus, the dynamic distance classifier computes decision boundaries (d) and their associated probabil- 
55 ities for correct classification (p) for each feature under consideration. The dynamic distance classifier relies 
on a predetermined tolerance zone T, for timely diagnosis of tool wear, and the fact that samples are 
sequential in time, to minimize the generation of false alarms. A false alarm is a diagnosis of a tool as being 
worn that is made while the tool is still good, i.e., before the tolerance zone. In step 6, by iteratively 
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updating the relevant parameters until the condition terminating the While loop has been satisfied, false 
alarms are avoided. Furthermore, p is weighted by s, the number of false alarms. 

Returning to Fig. 9 f the heuristic rules for computing the conditional probabilities when using the 
dynamic distance classifier differ from those described above with respect to tool life prediction and 

5 verification of tool life prediction. The decision boundaries generated by the dynamic classifier are assigned 
as the state boundaries for each feature node. Using the marginal probability of correct classification for 
each feature, some heuristic rules are used to automatically generate the overall knowledge base for the 
real-time expert system. In the influence diagram of Fig. 8, for the rms node 803, the decision boundaries 
for each category 0 to c, be denoted by ai ... a n , where n £ c+ 7, and the corresponding marginal 

70 probabilities for correct classification is denoted by r\ .. r n (n can be less than c + 7, because the decision 
boundaries of some the categories may be made to coincide if they lie close to each other. Similarly, by ... 
b n and r\ ... r n denote the decision boundaries and marginal probabilities for the Arms node 801. The node 
rms 803 is divided into n + 7 states and the node Arms is divided into m + 7 states. For example, node 803 
is in state 0 if x < ai , in state 1 if di £ x £ fe, ... in state n if x £ a n , where x is the rms feature value. The 

;s node 821 is in state 0 if the machining parameters correspond to category 0, ... , state c if the category is 
c. Also, the wear 7 node 817 is divided into two states, state 0 for a fresh tool and sfafe / for a worn tool. If 
k is the state of the category node 821, the critical state corresponding to the decision boundary and 
probability for the rms feature are k1 and r k1 , respectively, and for Arms, k2 and s k2 . 

Based on the above definitions, for each point (x x y) in the conditional cross-product space (Qsoz x 

20 Q8oi>) of nodes 803 and 801 , where x and y denote respectively the current state of nodes 803 and 801 , the 
following heuristic rules are used to compute the conditional probability that the wear7 node is in state 7, 
i.e., worn tool state. 

Conditional probability that node wear7 817 is in "worn drill" state for (x x y) € (flao3 x 0 8 oi^ 



25 



30 



35 



1.0 



kl 



k2 



k z + k 2 - x - y 



m ♦ n ♦ 2 
Jc t + k 2 - x - y 



m ♦ n «• 2 



x -y 



m * n ♦ 2 



if x 3t J^, y k 
if x i k lt y < k 2 

if x < k z , y * k 2 

if x < k lt y < k t 



If the computed quantities using the above equations lie beyond [0, 1], then the corresponding limit 
values are substituted for the conditional probabilities. Also, the above formulation can be extended to 
include more than three nodes as the predecessors of the result nodes like wear 7 817. 

40 Returning to Fig. 9, the dynamic distance classifier computes the decision boundaries and probabilities 
for correct classification. The latter are used in the heuristic rules described above to generate the 
conditional probabilities. These, combined with the decision boundaries and the relational structure, e.g. the 
influence diagram shown in Fig.8, are used to automatically generate the knowledge base. This step 
consists of writing the required information into the format for the knowledge base structure the Influence 

45 Diagram Expert System 154 expects. The generated knowledge base is used according to the method 
described above with respect to Figures 2 and 3. 

While this invention has been described with reference to illustrative embodiments, this description is 
not intended to be construed in a limiting sense. Various modifications and combinations of the illustrative 
embodiments, as well as other embodiments of the invention, will be apparent to persons skilled in the art 

50 upon reference to the description. It is therefore intended that the appended claims encompass any such 
modifications or embodiments. 

In particular and in a first aspect of the present inventions system for predicting the life of a cutting tool 
in an automated metal cutting machine, comprising: 



55 a plurality of sensors for measuring physical data relating to a cutting operation of said machine; and 

a computer connected to said sensors, comprising a feature extractor for transforming said physical 
data into feature values; 

a tool life predictor for arriving at a prediction of the life of said cutting tool, having 
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an influence diagram having an input level of feature nodes, intermediate nodes connected to said 
feature nodes and to a machining class node, an output node connected to said intermediate nodes; 
a first module for computing averages of feature values for each sensor; 

a second module responsive to said feature averages for making a classification of said feature average 
as indicative of a long life or a short life; and 

a third module connected to said influence diagram for combining said classifications to arrive at a 
classification and a conditional probability at each intermediate node and a classification at said output node 
which is the predicted life of said cutting tool. 

In a further aspect of the present invention, a system includes a tool wear diagnostic module responsive 
to said predicted life of said cutting tool and to said verification of said predicted life of said cutting tool, for 
arriving at a diagnosis of wear of said cutting tool, having 

an influence diagram having an input level of feature nodes, intermediate nodes connected to said 
feature nodes and to a machining class node, an output node connected to said intermediate nodes; 

a seventh module responsive to said feature values for making a classification of said feature values as 
indicative of a tool being worn or not worn; and 

a eighth module connected to said influence diagram for combining said classifications to arrive at a 
classification and a conditional probability at each intermediate node and a classification at said output node 
which is the diagnosis of wear of said cutting tool. 

20 Claims 

1- A system for predicting the life of a cutting tool in an automated metal cutting machine, comprising: 
a plurality of sensors for measuring physical data relating to a cutting operation of said machine; and 
a computer connected to said sensors, comprising a feature extractor for transforming said physical 
25 data into feature values; 

a tool life predictor responsive to said feature values to arrive at a predicted tool life of said cutting tool. 

2. The system of claim 1 and further comprising a tool life prediction verifier responsive to said predicted 
tool life and said feature values for arriving at a verification of said predicted tool life. 

30 

3. The system of claim 1 or claim 2 and further comprising a tool wear diagnostic module responsive to 
said predicted tool life, said verification and said feature values to continuously diagnose said tool for a 
worn out condition to arrive at a diagnosis of tool wear of said cutting tool. 

35 4. The system of any preceding claim wherein said tool life predictor, arrives at a prediction of the life of 
said cutting tool, comprises: 

a knowledge base having an influence diagram with an input level of feature nodes, intermediate nodes 
connected to said feature nodes and to a machining class node, and an output node connected to said 
intermediate nodes; 

40 a first module for computing averages of feature values for each sensor; 

a second module responsive to said feature averages for making a classification of said feature 
averages as indicative of a long life or a short life; and 

a third module connected to said influence diagram for combining said classifications to arrive at a 
classification and a conditional probability at each intermediate node and a classification at said output 
4$ node which is said predicted tool life of said cutting tool. 

5. The system of any of claims 2, 3 or 4 where said tool life prediction verifier comprises: 

a knowledge base having an influence diagram with an input level of feature nodes, intermediate nodes 
connected to said feature nodes and to a machining class node, and an output node connected to said 
50 intermediate nodes; 

a first module for computing maximum values of said feature values for each sensor; 

a second module responsive to said maximum values for making a classification of said maximum 

value as indicative of a long life or a short life; and 

a third module connected to said influence diagram for combining said classifications to arrive at a 
55 classification and a conditional probability at each intermediate node and a classification at said output 
node which is said verification of said predicted tool life of said cutting tool. 

6. The system of any of claims 3, 4 or 5 wherein said a tool wear diagnostic module comprises: 
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a knowledge base having an influence diagram with an input level of feature nodes, intermediate nodes 
connected to said feature nodes and to a machining class node, and an output node connected to said 
intermediate nodes; 

a first module responsive to said feature values for making a classification of said feature value as 
5 indicative of said tool being worn or not worn; and 

a second module connected to said influence diagram for combining said classifications to arrive at a 
classification and a conditional probability at each intermediate node and a classification at said output 
node which is said diagnosis of tool wear for said cutting tool. 

10 7. The system of any preceding claim, wherein said sensors comprise: 
a strain gage. 

8. The system of any preceding claim, wherein said sensors comprise: 
a spindle current sensor. 

;s 

9. The system of any of preceding claim, wherein said sensors comprise: 
a feed current sensor. 

10. The system of any preceding claim, wherein said sensors comprise: 
20 a dynamometer. 

11. The system of any preceding claim, wherein said sensors comprise: 
an accelerometer for sensing vibrations in said machine. 

25 12. A method of operating a metal cutting machine having a cutting tool, comprising the steps of: 

(a) attaching at least one sensor to said metal cutting machine; 

(b) commencing a metal cutting procedure; 

(c) observing wear of said cutting tool and noting when said cutting tool is worn out according to 
some criterion; 

30 (d) collecting data from said sensors; 

(e) selecting from said sensors features which influence the wear of said cutting tool; and 

(f) using said selected features for predicting the life of said tool. 

13. The method of operating a metal cutting machine of Claim 12, further comprising the step of: 
35 (g) verifying the predicted life at the mid-point of the predicted life. 

14. The method of operating a metal cutting machine of Claim 12 or 13, wherein step (e) further comprises 
the step of: 

(e.1) determining decision boundaries for the data collected in step (d). 

40 

15. The method of operating a metal cutting machine of Claim 14, wherein step (e) further comprises the 
steps of: 

(e.2) selecting machining parameters for categorizing the machining operation; and 
(e.3) building an influence diagram for modelling the wear of said cutting tool based on the data from 
45 said sensors. 

16. A method for learning in a real-time expert system, comprising the steps of: 

(a) creating a relational structure corresponding to the influences that a first and a second variable 
have on a third variable; 

so (b) obtaining data for said first variable, said second variable, and said third variable over some 

duration of time; 

(c) selecting a tolerance zone around a sample that divides classes; and 

(d) using a dynamic distance classifier, which is responsive to said tolerance zone, to determine 
decision boundaries for said first and second variables and probabilities associated with said 

55 decision boundaries, said probabilities indicating the confidence with which a decision boundary 

classifies values of said first and second variables with respect to the value of said third variable. 
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